New Value of m M /m e from Muonium Hyperfine Splitting 

Richard J. Hilfl 

Newman Laboratory of Nuclear Studies, Cornell University 
Ithaca, New York 14853 
(February 1, 2008) 

The complete contribution to the muonium hyperfine splitting of relative order a 3 (m e /m M ) In a 
is calculated. The result amounts to 0.013 kHZ, much smaller than suggested by a previous estimate, 
and leads to a 2a upward shift of the most precise value for the muon-electron mass ratio, with the 
error reduced by approximately 30%. Analogous contributions are calculated for the positronium 
hyperfine splitting: (217/90 - 171n2/3)m e (a 7 /7r) lna -1 « -0.32 MHz; the remaining theoretical 
uncertainty is well below experimental error, leaving discrepancies of 2.5a and 3.5a with the two 
most precise measurements. 

Precise measurement of the ground-state muonium (/i + e _ ) hyperfme-splitting (HFS), together with the correspond- 
ing theoretical analysis, provides a stringent test of bound state theory in Quantum Electrodynamics (QED), and 
allows a precise determination of the fundamental physical constants m M /m e and a. The most precise measurement 
> ! gives @: 

^ ; Ai/(Mu) eX pt. = 4 463 302.765(53) kHz [1.2 x 10~ 8 ]. (1) 

The corresponding theoretical prediction can be expressed as a series expansion in small parameters a w 1/137 and 
iTi e /m^ ks 1/207; terms involving logarithms, lna -1 rs \n.(m^/m e ) fss 5, also appear. At leading order in a, the 
, splitting is given by the Fermi energy (|]: 
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The complete splitting can be broken into the sum of terms § [|: 

A^(Mu)thoory = A^_D + A^ rad + A^ rcc + Al/ r _ r + A^ wca k + A^had- (3) 

■ Here D stands for Dirac, or relativistic corrections, while the other terms are from radiative, recoil, radiative-recoil, 
f> |j weak and hadronic contributions. 

JL i Currently, theory is limited by uncalculated or imprecisely-known terms in Af rec and Ai/ r _ r of order E F a 3 {m e /m^), 
some of which are enhanced by logarithmic factors; see Table I. This paper presents a calculation of terms of order 
r"| , EFCx 3 (m e /m IJ ,)\na, with results: 
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Numerically, these contributions give —0.034 + 0.047 = 0.013 kHz. Previous incomplete calculations |5| || sug- 
gested a contribution of —0.263(60) kHz. The main result of this paper is to show that in fact these contribu- 
tions are not as large as the previous estimates. Remaining theoretical uncertainty is dominated by terms of order 
E F (Za) 3 (m e /m^)ln(m^/m e ) (~ 0.06kHz); E F (Za) 3 (m e /m^) (~ 0.03kHz); and E F a{Za) 2 (m e /m^) (~ 0.03kHz). 
A discussion of the error due to still uncalculated terms is given at the end of the paper. 

Including the complete results of Eqs.([|),(|5|) does not significantly alter the theoretical prediction for the HFS 
in physical units, which is in good agreement with the experimental value, Eq.(|l]). Here we simply quote 
Ref. §(Eq.(D14)): 

A^(Mu)thcory = 4 463 302.67(27) kHz [6.1 x 10~ 8 ], (6) 

where the error is due mainly to the measured value of rnn/m e |Q. Likewise, the HFS determination of a is not 
significantly changed H. However, the new results in Eqs.(P,(@) represent a fractional shift of 6.2 x 10 -8 in the HFS, 



and hence also in the HFS determination of the mass ratio m M /m e . A recoil term of order E F (Za) 2 {m e /m ll ) 2 , which 
was not included in Az/ rec of Ref. contributes an additional 0.065(6) kHz or 1.5 x 10~ 8 [g. The mass ratio is then 
shifted from the value in Ref. ||(Eq.(161)) to become: 
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(m M /m e ) [Av(Mu)] = 206.768 2817(33) (24) (16) [2.1 x 1(T 8 ], (7) 

with the errors arising from uncertainty in A^thcory due to uncalculated terms, from Az/ expt ., and from the value of a, 
respectively. This represents a shift of 2.5a (in terms of the previous error 3.1 x 10~ 8 ), and a 30% reduction in error. 

The positronium (e + e~) HFS has also been measured precisely, though at present its interest is for testing our 
knowledge of QED bound states, as opposed to determining fundamental constants. The two most precise values are 
due to Mills and Bearman (Aj/(P) exp t. l, Ref- jf|) and Ritter et. al. (Ai/(P) oxpt . 2 , Ref. Jl(|): 

At/(Ps) cxp t. i = 203 387.5(1.6) MHz [7.9 x 10~ 6 ] (8) 

Ai/(Ps) oxpt . 2 = 203 389.10(74) MHz [3.6 x 10 -6 ]. (9) 

The theoretical expression is: 

Ai/(Ps)thoory = m e a 4 f C + Ci - + C 2 ia 2 In a" 1 + C 20 a 2 (10) 
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Including the known terms through C 2 o 11 yields A^(Ps) Q 2 = 203 392.93 MHz. Coefficient C32 = —7/8 has been 
known for some time fl2]| , and in this paper we calculate: 

C31 = 217/90- 17 In 2/3. (11) 

C32 and C31 contribute —0.91 MHZ and —0.32 MHz to Az/(Ps)theory, respectively, bringing the theoretical prediction 
to: 

A^(Ps)thcory = 203 391.70(20) MHz [1.0 x 10~ 6 ]. (12) 

The uncertainty of 0.20 MHz corresponds to a coefficient C 3 o ~ 4. For comparison, the numerical values of the other 
coefficients are: C = 0.58, d = -1.24, C 21 = 0.21, C 20 = -0.39, C32 = -0.88, C31 = -1.52. The discrepancy with 
experiment is significant: 2.5<r and 3.5cr for Eqs.(||) and (^|) respectively. As with the orthopositronium lifetime |l3| ] 
||14| , a true disagreement between experiment and the predictions of QED would have important consequences. 

The calculation is done in the framework of an effective quantum mechanical Hamiltonian theory 14 1, taking inputs 
from relativistic QED field theory and from (non-relativistic) NRQED field theory |l^] . The results to be derived for 
muonium can be translated directly to positronium by taking m M — > m ei and including the additional contributions 
from virtual e + e~ annihilation. 

The Hamiltonian can be decomposed into the sum: 

H = H + V 4 + V 5 + V 6 + V 7 , (13) 
where Ho is the unperturbed Hamiltonian for the Coulomb problem with reduced mass m r = m e m fJi /(m e + m^): 

f.-^_(£50. (14) 
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Potentials V4, V5, Vq and V7 give contributions to the energy of order ma 4 , ma 5 , etc. Since non-HFS operators will 
affect the HFS only in second- or higher-order perturbation theory, it follows that only the HFS parts of potentials Vq 
and Vj are necessary. Furthermore, any potential not contributing to S-states (in first or second order perturbation 
theory) may be neglected. 

We will write the potentials in terms of a list of standard operators: (q = I — k) 

(l\Oi\k) = \ (15) 

1 i ( P 2 Z a m r /2 i 

2 = 7r(Za)mf (2^ ~ ~ ~ J JL_^_£ P 

2m r r 

»|fe) = ^ffiM 



(i\o 4 \k) 



1 \l x k\ 2 
itf, q 2 



2 



(l\0 5 \k) = 


n(Za)-^— 


(l\O e \k) = 


4i„M 

ml m r 




1 p 


ir(Za) 


(l\O s \k) = 


1 \l x k\ 2 


m 2 . q 4 


(l\0 9 \k) = 


1 ( 

m r V 



3g ■ cr e g ■ 07 

<? 2 



Note that d = i 3 (r)/m^. 

Potential V4 is derived from tree- level NRQED diagrams with Fermi, Darwin, Kinetic, and Dipole vertices [jl6), and 
contains the leading relativistic corrections: 
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Renormalization constants cp ss = 1 + 0(a) are tabulated below. V5 gives the leading radiative corrections: 
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Using first order perturbation theory, potentials V4 and V5 correctly reproduce the complete S-state energy spectrum 
through 0(ma 5 ) 14 [l7[ |. The contribution from muon vacuum polarization is not relevant to our analysis, and has 
been excluded from V5. 

For Vq, only HFS terms are necessary. These again are taken directly from NRQED diagrams: 



V 6 = Air(Za) 
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Spin-Orbit, retardation, Time-Derivative, p'p, and Seagull interactions have been included. Additional local operator 
terms, of the form — V 2 5 3 (r) and {p 2 , 8 3 (r)} are not shown explicitly; these analytic terms do not generate factors of 
In a, and so are not relevant to the present analysis |Ts|| . 

The necessary renormalization constants have already been calculated [EBf JTol : 
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Here a e — a/2% + 0(a 2 ) is the electron anomalous magnetic moment. For c M , m M and Z 2 a are substituted for m e 
and a. 
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Potential Vj has no non-instantaneous HFS contribution coming from photon momenta q « ma 2 , a consequence of 
the fact that spin-dependent Ml multipole transitions vanish in the absence of relativistic effects, and are therefore 
suppressed. The remaining instantaneous part of Vj, from momenta q w ma, is fully determined by requiring that the 
Hamiltonian correctly reproduce the low-momentum expansion of the 1-loop photon-exchange scattering amplitude. 
Introducing a photon mass A, and ultraviolet cutoff A on photon momenta, the effective Hamiltonian (without Vj) 
gives pq|: 



2n(Za) 
3m e m u 
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where again analytic terms are not shown. The corresponding QED amplitude is (Fig. 1): 
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This result has been checked both in QED Feynman gauge, 
effective theory to match QED implies that 



and in NRQED Coulomb gauge pij|. Requiring the 
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Contributions to Vj having a dependence on m e , m^, a and Z different from Eq. (p2|) are ruled out by noticing that: 
(i) The non-recoil contributions are already present in V4, V5 and Vq (as we will soon verify), so that V7 contains no 
non-recoil piece; (ii) Masses can enter only as inverse powers l/m e and l/m M , and in particular not as l/(m e + m^). 
This latter result can be seen clearly using time ordered perturbation theory in NRQED: the NRQED vertices are all 
homogeneous in the masses, leaving only the energy denominators to consider; however, the energy denominators will 
all have the form l/(|q| + p 2 /m e + p\jm^), with photon momentum q and particle momenta p\, p2- (Contributions 
which are not simply iterations of lower-order potentials must have at least one photon in each intermediate state.) 
Such an expression, for q w p\ « P2 ~ ma, can be expanded in powers of p\jm e \q\, P2/m^\q\ — again homogeneous in 
the masses. Using (i) and (ii), the only possible parameter dependence which is symmetric in m e and m^ is that of 
Eq.©. 

Having completed the specification of the Hamiltonian, Eqs. (|T^) , ([li]) , jl^) , (|l7|) , (|T^) , (p2|) , we now use the usual 
expressions from Rayleigh-Schrodinger perturbation theory to solve for the energy shift: 



AE = (V 6 + V 7 ) + ( {Vi + V 5 )G{Vi + V 5 
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where G is the Coulomb Green's function with ground state pole removed, and (V) is the expectation value of V in 
the ground state of the unperturbed H$, Eq.(0). The logarithmic contributions of the necessary matrix elements are: 
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(25) 



(OoGOy) 10 (Za) x 
/X3r n V >— ln(Za) , (26) 

where the arrows signify that only logarithmic corrections, and in the case of (O9GO9), only the HFS part, are 
shown p2| . The pure recoil result for the HFS at order £'i?(2'a) 3 (m e /m AI ) contains the previously known \n 2 (Za) 
and ln(Za) ln(m M /m e ) contributions []l2| |l(| (^3); the new ln(Za) term is shown in Eq.(Q) (^4|. For radiative 
corrections at order Epa(Za) 2 , the non-recoil hi 2 (Za) and lu(Za) terms, and the recoil (m e /m^)hi 2 (Za) term [ p"2| , 
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agree with previous calculations. A part of the radiative-recoil single-logarithm corresponding to reduced mass factor 
m^./m e m il ~ (1 — 2m e /m tl ) was included previously Q 0; the complete contribution is given in Eq.(||). Numerical 
values are summarized in Table I. 

For positronium, there are additional interactions due to virtual annihilation of the electron and positron. The 
hard annihilation process is described by local operators, which by themselves cannot generate nonanalytic factors. 
So, for In a contributions, only second order perturbations involving V4 and V5 need be considered: 

SVa = f (| + ^) O x (27) 



a 
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8V4 gives the leading contribution from 1-photon annihilation. The first and second terms of 8V5 come from radia- 
tive corrections to 8V4, and from 2-photon virtual annihilation, respectively. O (ma In a) contributions from these 
annihilation operators are: 



A^ ann .(m e Q! 7 In a) = m e — In a 1 
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The non-annihilation contributions for positronium are obtained by taking the limit m M — > m e in the muonium 
analysis (making no expansion in m e /m M ); the combined result is given in Eq. (p"l|) ■ 

The previously most significant sources of error in the muonium HFS were Ai/ r _ r (0.104 kHz) and Ai/ rcc (0.060 kHz) 
H ; all other uncertainties are estimated below 0.010 kHz ^5|. By calculating the 0(EFa 3 (m e /m^) In a) contribution 
to Az/ r _ r , the uncertainty in this quantity should be reduced by a factor ~ lna -1 ~ 5; in fact, since there are still 
uncalculated terms at 0(Ep a 2 (Za)(TO e /m M ) ln(m M /m e )) |2(| and 0(EFa(Za) 2 m e /m^) we take this uncertainty 
as 0.040 kHz. The uncertainty in Av lec should remain approximately the same, since it is dominated by the still 
uncalculated terms of order 0(E F (Za) 3 (m e /m fJj )\n(m fJj /m e )) Q and 0(E F (Zaf (m e / m^)) pgQ . Thus we take 
0.070 kHz as an estimate of the total remaining theoretical error. 

In the final stages of the calculation, I received word from K. Melnikov and A. Yelkhovsky that they have also 
performed the calculation of a 3 In a terms, in a dimensional regularization approach p9]j . After a detailed comparison, 
we agree fully on the contributions in both muonium and positronium. The agreement found in different formalisms 
in two independent calculations lends strong support to the correctness of the results. 

This work was motivated in part by, and is an extension of, Ref. Jl4| . Many ideas used in the calculation originated 
with G. P. Lepage, who I thank for continued insights and encouragement during the present work. Thanks are also 
due to P. Labelle, and to K. Melnikov and A. Yelkhovsky for useful conversations. This work was supported by a 
grant from the National Science Foundation. 





Ref. § (kHz) 


present paper (kHz) 


(Zaf In* (Za) 
(Zaf In (Za) ln(m M /m e ) 
(Za) i \n(Za) 
(Zaf ln(m M /m e ) 
(Zaf 


-0.043 
-0.210 
-0.257(f) 

0.107(30) 


-0.034 

-0.035 (*) @ 


a(Zaf \n A (Za) 
a(Za) 2 In (Za) 
a(Zaf 
Z 2 a(Zaf In (Za) 


0.344 

-0.008 (*) 
-0.107(30) 


0.034 
0.013 


a i (Za) ln a (m M /m e ) 
a 2 (Za) ln 2 (m M /m e ) 
a 2 (Za) ln(m M /m e ) 
a 2 (Za) 


-0.055 
0.010 
0.009 (*) 





TABLE I. Contributions of order EFa 3 (m e /m ll ) to the muonium HFS. The second column lists the contributions used in 
Ref. [3]; the third column gives new or modified values from the present paper. Asterisks denote partial results. 
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Figure 1: Feynman diagrams for the QED amplitude necessary to determine V7. 



